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FORF.UOSD 


The results of an analytical Investigation of nonlinear, nonplanar analysis 
methods applicable to both thin-jet s»deling and thick-jet suxleling of STOL 
aircraft ground-effect phenomenon are presented. The study consisted of 
two concurrent tasks which are reported on separately as Parts I and II. 

Part I. Nonplanar. Honlinear Wing/Jet Lifting Surface Method . The objec- 
tive of this task was to extend the Douglas Nonplanar Lifting Systems 
program to include powered-lift wings having thin jets of varying strength 
for both part and full span arrangements and to analyze various configura- 
tions in ground effect. 

Part II. Nonplanar. Nonlinear Method Applicable to Three-Dimensional Jets 
of Finite Thickness . The objective of this task was to apply the NASA Ames 
Research Center Potential Flow Analysis to power-off ground effect cases 
and recommend procedures for developing a thick-jet analysis method. 

This study, conducted by the Technology Programs Section, Aerodynamics 
Subdivision of the Douglas Aircraft Company, was sponsored by the NASA 
Ames Research Center under Contract NAS2-9319. Dr. C. A. Shollenberger 
served as principal Investigator for the study under the technical direc- 
tion of Mr. D. N. Smyth. The NASA project engineer was Mr. David Koenig 
of the Large Scale Aerodynamics Branch. 

The contributions of Mr. M. I. Goldhammer, who served as the principal 
investigator on Part I during the early stages of the study, are greatly 
appreciated. His previous work on the development of the Nonplanar Lifting 
Systems program contributed significantly to the present work. The 
assistance of Mr. D. H. Neuhart in preparing the input and running many of 
the cases is also appreciated. 

The authors also gratefully acknowledge the contributions of Dr. R. T. Medan 
of the NASA Ames Research Center for his assistance and support in the 
application of the NASA Ames Research Center Potential Flow Analysis to 
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Th« influence of gi^ound proniaic^ on chree-dtaeiwlonel lift eyeteae hee 
been previouely etudied, ee reported in Kefereaee 1. ueint the Donilen 
Monplener Lifting Syeteae prograa. HefereMt 1 preaente calculated ateady 
and quaei-ateady aerodynaaic eharacteriatica of varioue «inga operatii^ la 
ground effect. Selected configurationa froa Reference 1 have been analysed 
in the present study using the NASA Aaes Research Center Potential Flow 
Analysis (POTFAN) in order to exaaine POTFAM's applicability to grotmd 
effects probleoM and to provide a second evaluation of the isqtortrace of 
ground proxiaity on some basic %rlng planforas. 

The consBonly addressed aircraft ground effect problem considers a lift sys- 
tem in motion at a constwit altitude above a ground plane. This steady- 
state representation, illustrated in Figure la. is an idealisation of the 
approach or departure of an aircraft frem a rumray. The ccmplete tmateady 
representation of aircraft motion near a ground plane, depicted in Figure 
Ic. includes the time history of the lift system which is manifested by 
varying strength shed trailing vorticity. Additionally, the unsteady model 
employs the proper geotaetric relationship between aircraft attitude. 0. 
aircraft angle of attack, a, and flight path angle, y. The quasi-steady 
representation of Figure lb Includes the geometrical relationships of the 
unsteady representation but ignores the lift system time history. Therefore, 
compared with the steady representation, the quasi-steady model provides an 
Improved description of the transient ground effect problem with only a 
minimal increase in computational ciXBplexlty. 

Results of two-dimensional vortex lattice calculations are presented in the 
next section to compare vortex lattice methodology, used in the present 
application of POTFAN, with exact results. Furthermore, the t%K>-dimenslonal 
trends for quasi-steady ground effect cases are presented for centrast with 
predictions for three-dimensional cases. The principal emphasis of the 
present study is examination of the three-dimensional POTFAN results for 
wings operating in ground effect as reported in Section 3. Finally, 
development of the capability within POTFAN to analyse propulsive influences 
on lift systems has been investigated and is discussed in Section 4. 
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2. TMO'DIMENSIOliAL VORTEX LATTICE C/O.CUUT10NS 


A slople airfoil analysis tMthod avploylng 41scrata vortieaa and ecmirantional 
vortex lattice eollocati(m rules was constructed to atqtplwMnt the three- 
disensional ground-effect calculations. Although a tvo-dlamsional analyaia 
option is available within the POTFAN eyat«s» the present two-disMaslonal 
results were obtained frou a rudis»ntary couputer prograu written specifically 
for this task. The two-diaensional study was initiated to investigate the 
suitability of vortex lattice oethodology used in the present three-dinensional 
POTFAN calculations to analyse wings oparating in ground effect. A secondary 
objective of the two-dimensional calculations was determination of aerodynamic 
characteristics for quasi-steady ascent and descent relative to the steady- 
state values. 

The two-dimensional analysis employed was a simple finite-elemnt formulation 
with equal length elements distributed along the airfoil chord. Each element 
was represented by a point vortex at a location one quarter of the element 
length aft from the forward element edge. The inviscid condition of tangential 
flow was imposed at an element control point which was located three-fourths 
of the element length from each element forward edge. All nonplanar and non- 
linear aspects of the airfoil geometry and velocity evaluation were Included 
In application of the airfoil boundary conditions and force evaluation. Actual 
image airfoils were employed in the analysis, rather than image or symmetry 
options, to avoid any possible errors or ambiguity in simulating the ground 
plane. Airfoil loadings were determined by calculating the force on each 
elementary vortex through application of the Kutta-Joukowski Law for the 
force on a point vortex. The flow velocities employed in this evaluation 
were the local velocity at each point vortex excluding the self-induced con- 
tribution. Formulation of this two-dimensional analysis is closely analogous 
to the three-dimensional POTFAN method as applied in Section 3 to analyze wings 
in ground effect. 

Calculated results obtained using the two-dimensional vortex lattice method 
are given in Figure 2 for steady and quasi-steady motion of a flat-plate air- 
foil near a ground plane. The quasi-steady cases are approximations to the 
complete unsteady ground approach or departure problem as described in the 


2 


prcvlotM sccticm. In Fi^rs 2 tha indicate tha Mtval caMt eon^utad 

during the study. Also, exact ualues of lift au^Mstatien froft Keferenee 2 
are shown in Figure 2 for the steady flight case. 

the lift aupumtation ratios predicted hy the vortex lattice aethod agree 
closely with the exact values of Reference 2 for the range of airfoil height 
above ground plane values considcMd. This agreMsnt provides confldei^c 
that vortex lattice nethodology is capable of predicting aarodynaaic grotmd 
effect. There are no exact values for craparison with lift augaentation 
ratios calculated for quaai-steady flight. However . Figure 10 of Reference 1 
indicates quasi-steady descent case exhibits greater lift augaentation than 
the steady case. This trend is in agreeaent with the present two-diaMnslonal 
calculations. It is apparent froa the present two-diaensional stisly that the 
quasi-steady ascent lift is generally predicted to be lower than the steady 
case for the flat place airfoil. 


3 


3. TntEE-DDOIISXOll/a. POTFAN CAZtCUUTIOIIS OF HZIIG8 » GiOUIlD 8FFECT 


This ssctlon prsssnts asrodpnnlc ehsrsctarlstlcs pradictsd by ths NASA Absj 
Rssasrch Ccntsr Potsntlsl Flow Analysis (FOTFAN) for various vii^s oparafiiim 
nasr a ground plana. All of ths prssantly eonsidarad erafigurations wara 
praviously analyssd using tha Douglas Nonplsnar Lifting Systans prograa as 
rsportsd in Rafsrsncs 1 and updatsd in Rafsranca 3. (All rasults curraatly 
presented ere from Reference 3.) Thus* the present results provide a coa* 
parison between the predictions of the two nathods and an additional evalu- 
ation of the significance of ground effect on lift systea characteristics. 

The POTFAM method is s potential flow singularity, finite element analysis 
that is Implemented by a series of modular coBq>uter programs which perform 
the various analysis functions. Tm of the program modules are described in 
References 4 and 5. Additional POTFAN documentation is provided by the exten- 
sive use of comments within the program code. The POTFAN program structure 
is exceptionally general in many respects including geoswtry input, singu- 
larity type and distribution, image options, and equation solving techniques. 

In the present study only a small number of POTFAN 's capabilities were employed. 

All POTFAN calculations discussed in the present section used vortex lattice 
collocation rules for vortex filament ai.d element control point locations. 

The potential singularity type selected for the POTFAN calculations was a 
lattice of discrete vortex filaments which can be equivalently formulated in 
terms of constant strength doublet distributions. The doublet formulation is 
employed within POTFAN and is a computationally efficient arrangement for 
application of vortex lattice rules. The element spacing over the wing sys- 
tems analyzed using POTFAN was generally uniform except for configurations 
where flap geometry required unequal element sizes. Also, a one-fourth elcMnt 
tip inset option was selected to improve solution convergence with respect to 
the nundier of %ri.ng elements. The undeformed trailing vortex wake position 
specified in the POTFAN input was always parallel to the ground plane which is 
not necessarily in the same direction as the freestream velocity for quasi- 
steady ascent or descent cases. 
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In orter to •••uro • ■uffieiont nunl»or of olMwnta ooro to ttMlyto Mch 
eimflguratlony « ••rloo of caleulatiaMi «ioro porfona^ for ««eh trtef with 
liier«Mii 4 fuaboro of ■aalyolo alnrato. The nuii>«r of oliawnto to bo «q»loyod 
MM Mloctod on the bants of smII variation In lift and drag praUetlona for 
incraaaiki nonbar of analysla alananta. Tba nudbar of apanwlaa at^ diordwlaa 
dlvtalona anployad for roiFAM conputatlona ara Indlcatad on aach figure shoir> 

Ing POTFAN results. 

The lift aupientation resulting frra steady notion near a groui^ plane la given 
in Figure 3 for a ractat^ular plMfom wing with aspect rutlo equal to 6.0. 

Several nethods were esq>loyed to analyse this configuration thereby providing 
an indication of the variation ancmg predicted ground effect Mrodynanlc char- 
acteristics. Hie lift predicted by the Dui^las Meunaiui Prograa (Reference 6) 
is taken from Figure 23 of Reference 1 and pertains to a 12-percent thick 
NACA 0012 airfoil section. All other results presented on Figure 3 are for sero 
thickness flat plate alrfjil sections. The Douglas Noi^lanar Lifting Systesa 
results are the revised values of Inference 3. POTFAM calculations were gener- 
ally conducted with ten spanwlse stui four chordwiae divisions but s(mm calcula- 
tions were perfomed wit’i a IS by 6 distribution to evaluate element spacing 
sensitivity. Finally, the Douglas Jet Wing Fuselage Program described in Refer- 
ence 7, was used to analyze the wing at two ground heights. The Douglas Ronplanar 
Lifting Systems, POTFAN, and Jet Wing Fuselage Program results for this aspect 
ratio 6.0 wing are all in close agreement and indicate smaller lift augmentation 
ratios than the Douglas Neumann Program results for all ground heights. In 
ground effect, tl^ influence of section thicteess (included in the Douglas 
Neumann Program calculationa) is comnonly assuMd to produce a negative lift 
or "such dotm" cosipared to zero thickness wing section lift. Hoimver, this 
thickness effect is not a rigorous result which cm be applied with confidence 
to all configurations. The Douglas Nonnlanar Lifting Systems, POTFAN and Jet 
Wing Fuselage Program results indicate an unexpected opposite trend from tiw 
"such down" concept. 

A second set of results for a rectangular wing operating in ground effect is 
presented in Figure 4 where lift curve slope is given as a function of leading 
edge height above the ground plane. The calculated lift slopes are for zero 
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thickness airfoil sections (coopered to the 22-percent thick experlaentel 
section) eml ere hesed on en esstnqption of s liiuier lift ilope between Mro 
end one degree engle of etteck. Ae in the cese of the espect retio 6.0 wing 
results of Figure 3* the Dougles Nonplener Lifting SystMM progm predicts 
slightly higher lift then tiie POOTAM veluee for this espect retio A.O wing. 

Also Figure 4 indicetes very close egre«sent between experiMntslly neesured 
lift slope and the POTFAM predictions. It should be noted that the CM^ined 
effects of wing thickness md fluid viscosity are m>t considered in either 
of the calculated results. 

The POTFAM results for a third rectangular planform wing are given in Figures 
S and 6. These figures correspond to the Douglaa Nonplanar Lifting Systeas 
results given by Figures 14 and IS of Reference 3. Figure 5 applies to a wing 
with flat plate airfoil section while Figure 6 corresporais to a zero thlckiMss 
airfoil section with a 40-percent chord full-span flap which is deflected 60 
degrees. Both figures include steady-state lift and drag valuea as well as 
quasi-steady ascent and descent values. The quasi-steady approximation to the 
fully dynatnlc transient ground effect problem employs the wing-ground plane 
geometry of the full dynamic case but ignores the lifting system w>tion time- 
history (see Figure 1). 

Comparison of the pre'.enc results with those of Reference 3 yields no sub- 
stantial differences in predicted aerodynamic characteristics. The quasi-steady 
analyses Indicate a greater lift augmentation than the steady-state case for 
descent and a smaller than the steady-state ratio for ascent. This trend in 
quasi-steady results are observed for both negative aix] positive lift incre- 
ments due to ground effects. This ordering of the ascent and descent values 
of lift relative to the steady-case results is consistent with the t%ro-dimensional 
vortex lattice results reported in Figure 2. 

Lift Induced drag predicted by POTFAN and the Nonplanar Lifting Systems program 
follows a regular trend and Indicates that the quasi-steady ascent drag ratios 
are consistently higher than the steady drag values while descent drag ratios 
are less than steady values. Application of vortex lattice methods, such as 
the present use of POTFAN. to calculate drag on wings with highly deflected 
flaps requires care in ensuring solution convergence and. typically, vortex 
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lactic* M*tho4* ov*r«*ciB*t* drag w^ltudaa CM^rad to ai^rUMiata nd 
other analytical iMthoda. The praaantly ca^loyad POTFAM force avaluatitm 
eetlw}d approxiJMtea the actual leading edge thruot tem by a dlatributlon 
of atreowla* fore* reaulting froe local application of tlM Kutta-JmikoraU 
Law for the force on a vortex filaaant. Although quite uMful, vertex 
lattice induced drag calculation* euat be viewed with caution. 

In order to axaaln* POTFAN ground effect calculation* for a «>ra cMiplax 
wing than the above exanplea, a typical tranaport wing «ma analyied with 
takeoff flap deflection. Figure 7, which correapMida to Figure 18 of Refer** 
en < 3, give* the incrciBent In angle of attack reaulting free notion near ^* 
ground aa a ftmction of wing lift coefficient. In thla caae the height of 
the wing above the ground ia apecified to be the height at idilch the aircraft 
main landing gear contact* the ground. Calculated POTFAM and Douglaa Nonplanar 
Lifting Syateaa reaulta are given in Figure 7 aa well aa unpubliahed vdnd 
tunnel data for thla tranaport configuration. POTFAN predlctima of Increoent 
in angle of attack are geiwrally lower than the value* eatiaated by the Douglaa 
Nonplanar Lifting Syatena program or the wind tunnel neaaurraent* except at 
low lift coefficient*. POTFAN calculation* with a larger nuiri>er of chordwiae 
wing eieaenta produced aubatantially the aame reaulta aa the 14 apanwlae and 
8 chordwiae dlviaiona uaed to obtain the reaulta of Figure 7. 
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4. DBVELOPHBIT OF THICK-JEI Aiya.T8X8 WITHXM POTFAM 


Currently tte ^4 Am* Keiearch Center Fotentiel Flow Anelysia (POTF^) la 
applicable only to alnple emircr Flow Mid tiMrefore it cmiot directly 
eatinate the nutual influence between aircraft propulaitm Mid lift ayatana. 
Several developMent optiems are availi^le to facilitate eatination of power 
effect! by P0TFAH. TlMae optiona include replacawnt of propuleive jeta by 
aolid boundariea, deflected jet repreMntationa, Jet flap tectaiiiquea and jet 
with finite thickncaa oodela. Each of tbeae optiona haa proven useful for 
analyaie of certain configuratiMis but often have linited generality. The 
thick-jet nodal propoaed in this aectiem ia a generally applicable jet repre- 
sentation within tiM finite eleaent MialysiSt potential flow femulatimi of 
POTFAH. 


Modeling of propulsive jets originated with propeller slipstrean analyses st^h 
as References 9, 10* and 11. Typically t)Mse aethods enploy linearised 
boundary conditiona at an assuaed jet location. Refinenent of propeller alip- 
streea aerodynaaic prediction sMthods involved ii^roved solution techniques 
and greater coo^lexlty of configurations idiich could be analysed within the 
fraaework of linearised theory. Achieveneot of STOL aircraft porforasnee 
through application of strong interaction between propulsion and lift systeas 
has stimulated new interest in jet interaction prediction beyond the scope of 
propeller slipstream analysis. Although models eaploying the jet flap ideali- 
sation of Reference 12 have proven useful for estiaating powered lift charac- 
teristics, the initial jet angle and aosmitua nuet be specified. Also, the 
jet flap model describes a thin high energy jet exhausting near the wing trail- 
ing edge and, therefore, is not directly appliceble to sons powered lift 
proposals such as upper surface blowing. 

In an effort to obtain a more useful jet interaction analysia, a thick-jet 
model is discussed below ’Jhich Is based on the analysis of References 7 and 
13. This Jet model is c>»patible trith the potential flow singularities avail- 
able within POTFAN and will provide a basis for discussion of factors which 
will require attention regardless of the jet aodel ultiaately selected. Since 
the purpose of the present study is to oaks general recowsMndations of 


procedurag for dc^loplng a thick-jct capability witbto FOTfAH^ etmaidara* 
tiona iritich are apeeific to a particular exiatiag MttM»d i^ll be avoided* 

Tbe tbick^Jet model auggeated preaently it aa invitcid* iacoapreaaible 
idealiaaticm of real jet flow. Correetiona for real fluid ef facta » auch aa 
Jot eotrainaent » cao be incorporated into tbe nodal aubaequently. Tba alter* 
native to tbe preaent fluid aaaui^tiona ia acceptance of a jet nodal wltb 
e9q>irically derived characterietica which reatrict tbe generality of tbe 
method. The area of intereat of nany thick-jet applicationa is the firat few 
jet dianefere downstrean of the jet origin. In this region viacous effects 
nay be secondary to the Jet ncNaentum Influence on lift systeas. Furthermore , 
jet models based on fundamental experircnts can lead to inconsistencies such 
as flow through solid surfaces. Therefore, the perfect fluid assumption will 
be employed, at least initially. 

A significant feature of finite element aerodynamic analyses, such as POTFAN, 
is the wide variety of configuration geometries that can be analysed. Since 
the proposed thick- jet analysis is also a finite element technique, few 
restrictions on the type of problem which can be analyzed should be necessary. 
For example, arbitrary initial jet shapes and locations are feasible. Also, 
jets may oilginate at free actuator disks (modeling propellers which do not 
impart swirl to the flow) or at duct trailing edges (forming an idealized jet 
engine). These basic jet capabilities combined with a finite element solid 
body analysis, such as POTFAN, can be employed to analyze very general con- 
figurations including externally blown flaps, deflected slipstreams and upper 
surface blown wings. 

The assumption of Invlscid fluid in all flow regions implies that, except for 
special areas (wing trailing vortex sheets, jet boundaries, etc.) the flow 
is irrotational and, therefore, a velocity potential can be defined. Addition- 
ally, the incompressible fluid assumption requires this velocity potential to 
satisfy Laplace's equation and consequently all of the usual potential flow 
solution techniques are available for this multi-energy jet interaction problem. 
For example, since the Laplace equation is linear, potential flow singularities 
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each Individually aatlsfylng Laplace's equatlon» may be auperia^iosed to obtain 
solutions for coiq>lex cases « 

Extension of POTFAN to multl*^nergy flow would not require alteration of tl» 
boundary conditions presently specified on solid bodies; however, application 
of the solid body boundary conditions would require inclusion of the Jet 
induced velocities to account for the Jet influence. The usual invlsc|d flow 
condition of tangential flow on solid surface boundaries would normally be 
ap,ilied for combined solid body and Jet computations; but other conditions, 
such as a specified flow normal to the body surface, could be accommodated as 
at present in POTFAN. 

Introduction of a thlck-Jet model into POTFAN will require boundary conditions 
to be specified on the Jet boundaries which separate regions of different 
energy. Invlscld jet fluid does not mix with the ambient flow and, conse- 
quently, one jet condition prohibits flow normal to the jet boundary surfaces. 
(In subsequent development of the jet model, an empirically specified flow 
normal to the jet boundary could be imposed to simulate entrainment of fluid 
by the jet.) 

Tn addition to the above jet kinematic boundary condition, a dynamic condition 
is necessary to insure that the jet boundary is unloaded. Since the jet 
boundary will not sustain a load, the static pressure is continuous across the 
boundary whereas the total pressure is discontinuous. These jet pressure con- 
ditions can be applied in niunerous forms, but it is apparent from the Bernoulli 
equation that the flow speed tangential to the jet boundary is discontinuous. 
This velocity behavior implies that the jet boundary is a vortex sheet and, as 
shown in Reference 13. the jet pressure condition requires that the cross- 
stream component of jet vortex sheet strength, y, to be given bv 



where &H is the total pressure difference between the two regions of flow, 
P is the fluid density and V is the mean local flow speed at the jet 
boundary. 
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A doublet sheet ie e^uivelent to e irortex sheet (see Reference for exs^ple) 
andt therefore* can also be employed to represent a Jet boundary, nieo* as 
shown in Reference 7* the Jet dynamic boundary condition c«a be eiqpressed by 





where u is the doublet strength and C ie tlw coordinate parallel to the 
local flow. The above two equations are Just two of the amny possible condi- 
tions on Jet singularity strengths which express the Jet dynamic condition. 

As discussed above* the thick-jet boundary is a surface with discontinuous 
tangential flow speed tdtlch is characteristic of a vortex sheet or equivalently 
a doublet sheet. Therefore* either vortex or doi^let distributions are 
appropriate potential singularities to represent the Jet boundaries within 
POTFAN. Continuous singularity distributions are desirable from an accuracy 
viewpoint but are likely to require prohibitively large computation time for 
a preliminary analysis effort. Also, since POTFAH (as applied in the present 
study) employs discrete singularity distributions to represent solid surfaces, 
a compatible singularity type specification for Jet boundaries is a discrete 
distribution. Reference 13 describes a Jet representation employing vortex 
filaments while Reference 7 discusses the constant strength doublet distri- 
bution Jet model illustrated in Figure 8. Either of these Jet boundary repre- 
sentations are compatible with the basic POTFAN singularity types and distri- 
butions. Other possible applicable Jet representations include ring or 
elliptical vortices such as used in Reference 14 or a linearized Jet model 
such as Reference 15. 

Within the proposed finite element analysis, the Jet boundary surface is 
divided into numerous elements. On each element an appropriate potential flow 
slngulari:;y Is placed to model the fluid behavior at the Jet boundary. Addi- 
tionally, a control point where the Jet kinematic and dynamic boundary condi- 
tions will be applied, is located on each element. Specific examples %fhich 
comply with these assumed characteristics are the Jet representations of 
References 7 and 13. 
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Probably tbe noat extenslva aodif Ication to POfFAN to facilltata 

iraltl*aaaxgy flow analysla involves developstent of a procedure for the appli- 
cation of the Jet and solid body boundary conditions to obtain a eonbined 
flovfield solution# The principal difficulty to be anticipated In the solu- 
tion process is introduced by the initially unknown Jet boundary location. 

To overcome this difficulty* an iterative solution scheme is required which 
successively approximates the Jet position while adjusting the strengths of 
the singularities representing the Jets and solid bodies. The final result 
of the solution process is a Jet boimdary location and set of Jet/body 
singularity strength values which simultaneously satisfy the boundary condi- 
tions prescribed above. Two iterative solution techniques will be discussed 
briefly. The first scheme divides the process of adjusting the Jet location 
and redetermining the singularity strengths into separate steps. In contrast* 
the second technique attempts to both move the jet boundary and reevaluate the 
singularity values simultaneously. 

The first iterative solution scheme, illustrated by the flow chart of Figure 9, 
is basically the procedure of Reference 13 and is initiated by an approxima- 
tion or "guess" for the jet position and singularity strengths. The matrix of 
aerodynamic Influence coefficients which express the solid body singularity 
distribution Influence of every body clen«nt on each body control point is 
then evaluated. This matrix, when multiplied by the vector of body singularity 
strengths, yields the total body Influence at each body control point. The 
influence coefficient matrix is next inverted (or equivalently trlangularized) 
so that the equations expressing the body boundary conditions can be repeatedly 
solved in an efficient manner. 

With these preliminary steps completed, the iterative process begins with cal- 
culation of the freestream and jet singularity influence at each body control 
point. The freestream and jet Influence are the nonhomogeneous part of the 
equations expressing the body boundary conditions. The body singularity 
strengths are then determined by multiplying the Inverted body influence 
coefficients matrix by the nonhomogeneous terms of the solid body boundary 
condition equation. 
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The newly detenained body singularity etrengtha penalt the flow velocities at 
each Jet boundary eleiaent control point to be detenained. Then the Jot is 
repositioned by forcing each eleaent to be tangent to local flowdet kinesatic 
condition) and the Jet singularity strengths are re-evaluated to satisfy the 
Jet pressure requireiaent (Jet dynaaic boundary condition). At this point in 
the solution cycle, one iteration has been coagtleted and a convergence test to 
detenalne the change in Jet/body properties is perfonaed. If a significant 
change in these properties has occurred during the s»st recent cycle, the steps 
beginning with calculation of Jet and freestreas influence on wing control 
points are repeated. A converged solution, very closely satisfying all Jet and 
body boundary conditions, is obtained tdien the variatic.. of Jet and body prop- 
erties bet%feen successive iterative cycles becoi^s saiall. At this occurrence, 
the iterative process is terminated and t\e pressures, forces, and other desired 
flowfield qualities are calculated. 

The second iterative solution procedure to be discussed is an adaptation of the 
method employed in Reference 16 to analyze wings with attached free vortex 
sheets. Bristow in Reference 17 applied a similar technique to calculate air- 
foil shapes with prescribed pressure distributions. Figure 10 gives a flow 
chart of this second procedure. Once again an initial starting solution or 
approximation is required to initialize the solution procedure. This starting 
solution must Include the Jet location as well as singularity strengths of 
each jet and body element. Next the Jacobian matrix of the equations expres- 
sing the solid body ami jet boundary conditions is calculated. This Jacobian 
matrix wl. a multiplied by the perturbations in the solution variables (singu- 
larity strengths and jet coordinates) is equal to the error betmen the present 
solution and the exact compliance with the boundary conditions. 

The error between the desired and present solution is then evaluated permitting 
the perturbations in solution variables to be determined. This process deter- 
mines updates to all singularity strengths and jet boundary coordinates in a 
single step. Next, these updates are employed to reposition the Jet boundaries 
and evaluate the jet/body singularity strengths. A convergence test checks the 
progress of the solution to determine if additional iterative cycles are 
required. Finally, desired flowfield properties are calculated from the con- 
verged solution. 
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Either of these two solution schemes is suitable for Incorporation into POTF^ 
for application of jet and body boundary conditions to obtain nulti-energy fl<^ 
field solutions. The first sch«ae offers the following advantages co^iared to 
the second method: 

1. Smaller matrices are involved since only the body boundary conditions 
are expressed in matrix equation fora. 

2. An efficient "inverse" matrix solution technique can be employed since 
the matrix of influence coefficients is only calculated once. 

3. Most required calculated quantities (influence coefficients* flow 
velocities* etc.) are already evaluated within POTFAN. 

The second iterative schene possibly provides benefits In rate of convergence 
and assurance of obtaining converged solutions. The potential benefits of the 
second schenm appear most certain when the solution scheme is employed in con** 
junction with higher-order singularity distributions (such as employed in 
References 16 and 17) rather than with the relatively crude distributions 
currently used in most POTFAN applications. The low-order of POTFAN singular- 
ity distributions combined with the three advantages listed above suggest 
that the first solution technique is probably most suitable for jet analysis 
within POTFAN. 

At the completion of the iterative solution process* the calculated results 
include the jet boundary positions and the strengths of singularities repre- 
senting bodies and jets. These quantities are usually of secondary interest 
compared to the body pressures and flowfleld velocities. However, the calcu- 
lated jet position is useful for determining the realism of calculated solu- 
tions and indicating problem areas for solutions that fail to converge. 

Graphic presentation of the calculated jet location is a valuable aid in 
assessing solutions and would be a logical addition to the existing POTFAN 
graphics capability. 

The flow velocity at any point in the jet/body flowfleld can be evaluated from 
converged solutions by summaticn of the jet and body singularity induction 
with the freestream influence. Consequently* the pressures on bodies can be 
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detemiiied using the local i^loclty on the body surface including tlM jat 
influence. Presancly» one POTFAM option determines the force on body olaunits 
by repeatedly ^plying the Kutta-Jotdtotraki Iisv for the force on a vortex fils- 
aent. In this procedure the velocity at each vortax filanent midpoint is 
oiployed and nusMrical difficulties cou&d be encountered if the filaamt mid- 
point is in close proximity to a Jet boundary which Is representMl by discrete 
singularities. In Reference 7 this difficulty is avoided by using the velocity 
only at element aldpolnts which can usually be selected to avoid close 
encounters with concentrated jet singularities. Some care may be required 
within POTFAN to precli^e numerical problems in the force evaluation procedure. 

The above discussion described general properties of a thick-jet model with 
only minor eaphasls on specific considerations relative to POTFAN. The rmalnder 

of this section will examine areas of POTFAN which %rlll require n»dification 

I 

to provide a thick-jet analysis capability. 


Provisions will be necessary within POTFAN to describe the jet characteristics 
and the jet properties required to initialise the iterative solution scheme. 

The existing POTFAN input program is very general and can be simply adapted to 
include jet related input. Jet properties which must be input Include the Jet 
strength (probably characterized by the total pressure difference relative to 
ambient flow conditions), jet origin coordinates and whether the jet originates 
at a free actuator disk or at the trailing edge of a duct. Additionally, the 
finite element properties such as the number of circumferential and stream- 
wise elements, the element spacing and the do«mstream extent of the jet calcu- 
lation before a downstream representation is employed, must be specified. 

Also, the starting or "guess" values for jet singularity strengths and jet 
boundary locations must be approximated to Initiate the solution process. 
Prevision for optional use of linearized or simplified solutions to approximate 

^ the jet singularity strengths is desirable since these strengths are not easily 

I i 

I determined Intuitively by users. Ease in approximating the initial jet posi- 

I 

P { tlon Is an Important consideration in user acceptance of the miltl-*energy flow 

h:T~~ ^ 

If i analysis. Finally, some input will be required to specify the ftmetions to be 

pet formed and techniques to be applied during the iterative solution process. 
For example, the number of iterations to be performed, convergence criteria, 
relaxation of solution variables, and sublterative cycles must be specified. 



Also, input provisions snist be stsde for optionel functions such ss force 
evalustion, offbody velocity evaluation, graphic presentation of solutions 
or restart of the iterative process fr<m a previous, partially converged 
solution. 

Prior to entering the iterati^ solution cycle sosw prelisetnar:; calcula- 
tions are required. The extent of these calculations depends on the iter- 
ative solution nethod selected. For exa^>le, for the schoe of Figure 9, 
the %ring influence coefficient natrix is evaliiated «id provisions ara SMde 
to solve the natrix equation enforcing the wing boundary conditions in an 
efficient nanner (inversion or triangularlsation) . Siallarly, in the nethod 
of Figure 10 elenents of the Jacobian natrix which are not altered by jet 
boundary position changes are calculated and stored. Many of these prelin- 
inary functions would be performed by the existing POTFAN influence coef- 
ficient evaluation and matrix solution programs with minor modification. 

Upon completion of the preliminary steps discussed above, the main iterative 
process is entered. The functions performed during this process are highly 
dependent on the iterative scheme selected. Basically, steps must be pro- 
vided to update body/jet singularity strengths and reposition the Jet 
boundaries. The jet reposition perhaps introduces the greatest difficulty 
in modifying POTFAN for nultl-energy flow analysis. Within the iterative 
scheme of Figure 9, the jet is repositioned using the jet kinematic boundary 
condition to align the jet surface with the local flow. In contrast, the 
method of Figure 10 requires the calculation of the elements of the Jacobian 
matrix which express the jet boundary conditions. In this case, the jet 
coordinates are part of the matrix equation solution. 

Most functions within the Iterative cycle, other than the jet reposition, are 
already available within POTFAN. For example, matrix solution, solid body 
boundary condition application, and velocity evaluation techniques are pre- 
sently operational. However, some modification to existing program code will 
be required, such as inclusion of jet induction in the velocity calculations. 

A necessary addition to the POTFAN structure for Toulti-energy flow analysis is 
an executive program to sunitor and control the iterative solution process. 


16 


The executive prograe direct! t)M fmctlona requested by the user Input end 
applies the numerical procedures specified. Integral to the executive progran 
Is calcoleclon of appropriate convergence test paraaeters on idilch execution 
decisions are based. The convergence parwBeters reflect the rate of change of 
the solution between Iterative steps as well as Indicate %ffaen a solution with 
acceptable coiq>llance with jet and body boundary conditions has be«i achieved. 
Also the executive progran could apply relaxation of the solution variables 
(jet singularity strengths» body slngularlt: strengthSt jet boundary coord* 
Inates) to Increase the rate of solution convergence or assist In obtaining 
convergence for difficult cases. 

Since It Is desirable to nalntaln POTFAM's present Independence of nachlne 
type, the executive progran should be written In FORTRAN rather than relying 
on manipulation of a particular machine control language. The modular POTFAN 
structure, with storage of the calculated result of each nodule in data sets, 
lends Itsrlf to control by an executive progran. 

The above discussion Identifies factors which must be considered to provide a 
thlck-Jet analysis capability td.thin POTFAN. This provides a basis to suggest 
a series of steps to undertake the modification to POTFAN. First, exanple 
computations should be completed using existing analysis codes with jet nodels, 
such as References 13, 14, and 15. These calculations would afford an exami- 
nation of the potential usefulness and possible difficulties of the candidate 
jet models and solution techniques. Next, the analysis components (such as 
Input modifications, jet reposition algorithm, and application of jet dynamic 
boundary condition) could be individually constructed and tested. Then It 
would be appropriate to Incorporate Into POTFAN the various analysis components 
along with an executive control program to direct the Iterative solution 
process. Finally, the resultant method would require validation and doctnsnta- 
tlon of the thick-jet analysis. 
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S. C0!iCLUDIIi6 REHARKS 


S«v«r«l conclttsioM «r« «ppar«it froa tiM pr«s«ic study on tlin appllc^lllty 

of ths Ames Research Center Potential Flow Analyais (POTFAB) to predict grotmd 

Influci^e and propulsive lift aysten aerodynaalc characteristics. The aost 

significant results are: 

1. Calculations for a tuo-dlaenalonal flat plate airfoil Indicate that 
vortex lattice analysis is useful for evaluation of ground InfliMnce 
on aerodynaalc characteristics. 

2. Two-dlaenslonal calculations employing the quasl**steady approximation 
indicate that a descending flat plate airfoil has greater lift augasnta- 
tlon. and an ascending airfoil has less lift augmentation, than the 
steady case. 

3. POTFAN calculations performed In the present study are in close agree- 
ment with the Douglas Nonplanar Lifting Systems results. Howsver, 

POTFAN generally indicated slightly lower lift augmentation, due to 
ground proximity, than the Douglas Nonplanar Lifting Systems method. 

4. The three-dimensional quasi-steady lift predictions agree with the 
ordering of ascent, steady flight and descent values found for the 
two-dimensional calculations. 

5. Fr(xn the limited number of cases examined, the POTFAN and Douglas Nonplanar 
Lifting Systems predictions generally appear to agree equally well with 
experimentally determined ground effect aerodynamic trends. However, 
effects of wing section thickness and fluid viscosity Introduce uncertainty 
in comparisons of theoretical and experimental ground effect results. 

6. POTFAN is suitable for extension for imiltl-energy flow analysis. Required 
modifications have been identified and systematic development steps have 
been suggested to provide a capability within POTFAN to analyse inter- 
actions between propulsion and lift systems. 
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Figure 3. Predicted Lift Variation with Ground Proximity 
of Rectangular Aspect Ratio 6.0 Wing 
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Figure 4. Experimental and Calculated Lift Slope 
Variation with Ground Height for 
Rectangular Wing. 
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Figure 5. Lift and Induced Drag of a Rectangular Wing 

in Ground Effect Predicted by POTFAN for Steady 
and Quasi-Steady Flight. 
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Figure 




RECTANGULAR WING 

FULL SPAN FUP 

FLAP DEFLECTION - 60 DE6S. 

FLAP CHORD » 0.4 TOTAL CHORD 

ASPECT RATIO -7.0 

ANGLE OF ATTACK - 10.0 DEGS. 

HEIGHT ABOVE GROUND MEASURED FR(W WING LEADING EDGE 
— O — POTFAN 15x7 STEADY Ct. » 3.6053 
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6. Lift and Induced Drag of a Rectangular Wing with 40-Percent 
Chord Flap in Ground Effect Predicted by POTFAN for Steady 
and Quasi -Steady Flight. 
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o WIND TUNNEL DATA 


NONPLANAR LIFTING SYSTQf PI^RAM 


POTFAN (14x8 eleaent-spacing) 


ASPECT RATIO -8.7 

TAPER RATIO - 0.2 

QUARTER CHORD SWEEP - 24 DECS. 

FLAP DEFLECTION - 5.0 DECS. Q 
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Figure 8. Jet Boundary Doublet Representation 






Figure 9. Jet Interectlon Anelyeie Iterative Solution Scheme Nuid>er 1. 
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Figure 10. Jer. Interaction Analyele Iterative Solution Scheme Nuaber 2 
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